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ABSTRACT 

Debris disks around main-sequence stars are believed to derive from planetesimal populations that 
have accreted at early epochs and survived possible planet formation processes. While debris disks 
must contain solids in a broad range of sizes — from big planetesimals down to tiny dust grains — 
debris disk observations are only sensitive to the dust end of the size distribution. CoUisional models 
of debris disks are needed to "climb up" the ladder of the coUisional cascade, from dust towards 
parent bodies, representing the main mass reservoir of the disks. We have used our coUisional code to 
generate five disks around a sun-like star, assuming planetesimal belts at 3, 10, 30, 100, and 200 AU 
with 10 times the Edgeworth-Kuiper-belt mass density, and to evolve them for 10 Gyr. Along with an 
appropriate scaling rule, this effectively yields a three-parametric set of reference disks (initial mass, 
location of planetesimal belt, age). For all the disks, we have generated spectral energy distributions 
(SEDs), assuming homogeneous spherical astrosilicatc dust grains. A comparison between generated 
and actually observed SEDs yields estimates of planetesimal properties (location, total mass etc.). As 
a test and a first application of this approach, we have selected five disks around sun-like stars with 
well-known SEDs. In four cases, we have reproduced the data with a linear combination of two disks 
from the grid (an "asteroid belt" at 3AU and an outer "Kuiper belt"); in one case a single, outer 
component was sufficient. The outer components are compatible with "large Kuiper belts" of 0.2-50 
earth masses (in the bodies up to 100 km in size) with radii of 100-200 AU. 

Subject headings: circumstellar matter — planetary systems: formation — Kuiper belt — stars: 
individual (HD 377, HD 70573, HD 72905, HD 107146, HD 141943) 



1. INTRODUCTION 

Since the IRAS disco very of the exc e ss infr ared emis- 
sion around Vega by lAumann et al.l ()1984l ). infrared 
surveys with IRAS, ISO, Spitzer, and other space- 
based and ground-based telescopes have shown the 
Vega pheno menon to b e common for main-se q uence 
stars (e.g. 'Mover et all l2004l: iBeic hman et aP 120051: 
Naiita, fc W iUiams 2005[Riekeet al.|l2005; Brvdc n et al l 
2p06; 'Sieglcr et al. 2007: 'Su et al." '2006; TrininK^Tal] 
2007 : Hillenbrand et al. 2008 : Trilling et al. 2008) . The 
observed excesses are attributed to circumstellar disks 
of second-generation dust, sustained by numerous plan- 
etesimals in orbit around the stars. Jostling collisions 
between planetesimals grind them all the way down to 
smallest dust grains which are then blown away by stellar 
radiation. While the bulk of such a debris disk's mass is 
hidden in invisible parent bodies, the observed luminosity 
is dominated by small particles at dust sizes. Hence the 
studies of dust emission have a potential to shed light 
onto the properties of parent planetesimal populations 
as well as planets that may shape them and, ultimately, 
onto the evolutionary history of circumstellar planetary 
systems. 

However, there is no direct way to infer the properties 
of invisible planetesimal populations from the observed 
dust emission. Dust and planetesimals can only be linked 
through models. First, dynamical models can be used to 
predict, for a given planetesimal family (mass, location, 
age, etc.), the distribution of dust. Such model s have 
become available in recen t years (e.g. [Thcbault et al.l 
120031 : iKrivov et all 120061 : iThebault fc Augereaul boOTt 



IWvatt et al.l 120071 : iLdime et all 12008"). After that, stan- 
dard thermal emission models will describe the resulting 
dust emission. Comparison of that emission to the one 
actually observed would then reveal the probable proper- 
ties of underlying, dust-producing planetesimal families. 

In this paper, we follow this approach and generate a 
set of hypothetical debris disks around G2 dwarfs with 
different ages (10 Myr - 10 Gyr), assuming debris dust to 
stem from planetesimal belts with different initial masses 
at different distances from the central star. For every set 
of these parameters, we simulate steady -state dust dis- 
tribu t ions with our coU isional code (Kri vov et al.l [20051 . 
120061: lL5hne et all [20081 ). This is different from a tradi- 
tional, "empirical" approach, in which dust distributions 
are postulated, usually in form of power laws, parame- 
terized by ranges and exponents that play a role of fit- 
ting parameters (e.g. Wolf & Hillenbrand 2003). Inter- 
estingly, replacing formal dust distributions with those 
coming out of dynamical modeling does not increase the 
number of fitting parameters. Just the opposite: the 
number of parameters reduces and those parameters that 
we keep free all have clear astrophysical meaning. The 
most important are lo cation of a parent p lanetesimal belt 
and its current mass (jWvatt et al.ll200'n) . 

Having produced a set of model debris disks, we com- 
pute thermal emission fiuxes in a wide range of wave- 
lengths from mid-infrared to millimeter. In so doing, we 
completely abandon simple blackbody or modified black- 
body calculations and solve a thermal balance equation 
instead. At this stage, we assume cor apact spherical 
grains composed of astronomical silicate ([Laor fc Draind 
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Il993f) and employ standard Mie calculations to compute 
dust opacities. Although this is still a noticeable simpli- 
fication, it represents a natural step towards considering 
realistic materials and using more involved methods of 
light scattering theory that we leave for subsequent pa- 
pers. 

As a test and a first application of the results, we re- 
interprete available observational data on a selection of 
disks around sun-like stars with well known spectral en- 
ergy distributions (SEDs). 

This paper is organized as follows. Section 2 describes 
the dynamical and thermal emission models. In section 3, 
a set of reference disks is introduced and the model pa- 
rameters are specified. Section 4 presents the modeling 
results for this set of disks: size and spatial distribution 
of dust, dust temperatures, and the generated SEDs. Ap- 
plication to selected observed disks is made in section 5. 
Section 6 summarizes the paper. 

2. MODEL 

2.1. Dynamical model 

To simulate the dust production by the planetesi- 
mal belt and the dynamical evolution of a disk, we use 
our coUisional code (ACE, Analysis of CoUisional Evo- 
lution). The code numerically solves the Boltzmann- 
Smoluchowski kinetic equation to evolve a disk of solids 
in a broad range of sizes (from smallest dust grains to 
planetesimals) , orbiting a primary in nearly Keplerian 
orbits (gravity -I- direct radiation pressure -I- drag forces) 
and experiencing disruptive and erosive (cratering) col- 
lisions. Collision outcomes are simulated with available 
material- and size-dependent scaling laws for fragmenta- 
tion and dispersal in both strength and gravity regime. 
The current version implements a 3-dimensional kinetic 
model, with masses, semimajor axes, and eccentricities 
as phase space variables. This approach automatically 
enables a study of the simultaneous evolution of mass, 
spatial, and velocity distribution of particles. The code 
is fast enough to easily follow the evolution of a debris 
disk over Gyr timescales. A detailed description of our 
approach, its numerical implementation, and astrophys- 
ical applications can be found in our previous p apers 
(iKrivov et al.ll200d[2005L[200llL5hne et al.l[200l . 
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For spherical dust grains with radius s and tempera- 
ture Tg we can calculate their distance r to the star under 
the assumption of thermal equilibrium as 




= dAQf^(s)i?A(rg) ■ 



(1) 



Here, i?* denotes the radius and F\^^,{T^,) the flux of the 
star with an effective temperature T* and Bx{Tg) the 
Planck function. The absorption efficiency Qf^^{s) is a 
function of wavelength A and particle size. 

We now consider a rotationally symmetric dust disk 
at a distance D from the observer. Denote by N{r, s) 
the surface number density of grains with radius s at a 
distance r from the star, so that N{r, s)ds is the number 
of grains with radii [s, s + ds] in a narrow annulus of 
radius r, divided by the surface area of that annulus. 
Then the specific flux emitted from the entire disk at a 



given wavelength can be calculated as 



^Hisk- / dr / dsi^A,disk(r,s) 
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xNir,s)Qris)BxiTg). 

3. REFERENCE DISKS 

3.1. Central star 

The parameters of the central star (mass and photo- 
spheric spectrum) affect both the dynamics of solids (by 
setting the scale of orbital velocities and determining the 
radiation pressure strength) and their thermal emission 
(by setting the dust grain temperatures). We take the 
Sun (a G2V dwarf with a solar metallicity) as a central 
star and calculate its p hotospheric spectrum w ith the 
NextGen grid of models (|Hauschildt et al.lll999( l. 

3.2. Forces 

In the dynamical model, we include central star's grav- 
ity and direct radiation pressure. We switch off the 
drag forces (both the Poynting-Robertson and stellar 
wind drag), which are of little importance for the optical 
dept hs in the range fro m ~ 10~^ to ~ lO ""^) cons idered 
here (|Artvmowic2jll997l : lKrivov et al.l[2000LlWvattl[2005h . 

3.3. Collisions 

The radii of solids in every modeled disk cover the in- 
terval from 0.1 /im to 100 km. The upper limit of 100 km 
is justified by the fact that planetesimal accretion models 
predict larger objects to have a steeper size distribution 
and thus to co ntribute less to the ma ss budget of a de- 
bris disk (e.g. 'Ken von fc Luuiri999b( l. To describe the 
coUisional outcomes, we make the same assumptions as 
in Lo hne et al.l (|2008r ). This applies, in particular, to 
the critical energy for disruption and dispersal, Q*u{s), 
as well as to the size distribution of fragme nts of an in- 
dividual collision. However, in contrast to lL5hne et al.l 
(2008) where only catastrophic collisions were taken into 
account, we include here cratering collisions as well. This 
is necessary, as cratering collisions alter the size distri- 
bution of d ust in the disk markedly, which shows up in 
the S EDs (jThebault et al.l l2003t iThebault fc Augereaul 
|2007( ). The actual model of cratering colli s ions u sed here 
is close to that bv IThebault fc Augereaul (2007). An es- 
sential difference is our assumption of a single power law 
for the size distribution of the fragments of an individ- 
ual collision instea d of the broken power law proposed 
originally in Theb anjlt et al.l (|2003f ). However, this dif- 
ference has little effect on the resulting size distribution 
in coUisional equilibrium. 

3.4. Optical properties of dust 

An important issue is a choice of grain composition 
and morphology. These affect both the dynamical model 
(through radiation pressure efficiency as well as bulk den- 
sity) and thermal emission model (through absorption 
efficiency). Here we assume compact spherical grains 
compos ed of astronomical sili cate (a.k.a. astrosilicate or 
astrosil. iLaor fc Draindll993[ ). similar to the MgFeSi04 
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olivine, wit h density of 3.3 g cm~'^ . Taking optical con- 
stants from iLaor fc Draind ()1993l ). we calculated radi- 
ation pressure efficiency Qpr and absor ption efficiency 
Qabs with a standard Mie routine (Boh ren fc HuffmanI 

To characterize the radiation pressure strength, it is 
custor nary to use the ra diation pressure to gravity ra- 
tio j3 ([Burns et al.lll979l ). which is independent of dis- 
tance from the star and, for a given star, only depends 
on Qpr and particle size. If grains that are small enough 
to respond to radiation pressure derive from collisions of 
larger objects in nearly circular orbits, they will get in or- 
bits with eccentricities e ~ /3/(l — /9). This implies that 
grains with f3 < 0.5 remain orbiting the star, whereas 
those with P > 0.5 leave the system in hyperbolic orbits. 
The (3 ratio for compact astrosil grains, computed from 
Qpr, is shown in Fig. [TJ The blowout hmit, /3 — 0.5, 
corresponds to the grain radius of s = 0.4 yum. Note that 
the tiniest astrosil grains (< 0.1 /im) would have (3 < 0.5 
again and thus could orbit the star in bound orbits. How- 
ever, the dynamics of these small motes would be sub- 
ject to a variety of effects (e.g. the Lorentz force) not 
included in our model, and their lifetimes may be short- 
ened by erosion processes (e.g. stellar wind sputtering). 
Altogether, we expect them to make little contribution 
to the thermal emission in the mid-IR to sub-mm. By 
setting the minimum radius of grains to 0.1 /im, we there- 
fore do not take into account these grains here. 
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Fig. 1. — Radiation pressure to gravity ratio /3 for astrosilicate 
grains as a function of their size. Horizontal lines at /3 = 0.5 and 
/3 = 1.0 show, which particles typically move in bound elliptic 
orbits, in hyperbolas, as well as in anomalous hyperbolas (open 
outward from the star). 



The spectral dependence of the absorption efficiency 
Qabs of different-sized astrosil spheres is depicted in 
Fig.H 

3.5. Parent planetesimal belts 

To have a representative set of "reference" debris disks 
around sun-like stars, we consider possible planetesimal 
rings centered at the semimajor axes of a = 3, 10, 30, 
100, and 200 AU from the primary. All five rings are as- 
sumed to have the same relative width initially (again, 
in terms of semimajor axis) of Aa/a — 0.2 (±0.1) and 
share the same semi-opening angle (the same as the max- 
imum orbital inclination of the objects) of e = 0.1 rad. 
The orbital eccentricities of planetesimals are then dis- 
tributed uniformly between 0.0 and 0.2, in accordance 
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Fig. 2. — Top: absorption efficiency of astrosilicate compact 
spherical grains as a function of wavelength for different grain sizes. 
Bottom: the spectrum of a G2 V star and the Planck curves for 
150 and 20 K (in arbitrary vertical scale) to indicate the spectral 
ranges most important for absorption and emission. 



with the standard equipartition condition. The initial 
(differential) mass distribution of all solids is given by a 
power law with the index 1.87, a value that accounts for 
the modification of the classical Dohnanyi's (|1969f ) 1.833 
thro ugh the size dependence of material strength (see, 
e.g.. lDurda fc Dermottlll997l ). 

The initial disk mass is taken to be IM0 (earth mass) 
for a 30 AU ring, roughly corresponding to ten (or slightly 
mor e) times the Edgeworth-Kuiper belt (EKB) mass 
(e.g. IGladman et al.l[200H IHahn fcMalhotral [20051) . For 
other parent ring locations, the initial mass is taken in 
such a way as to provide approximately the same spa- 
tial density of material. Since the circumference of a 
ring 27ra, its absolute width Aa, and vertical thickness 
2ae are all proportional to a, the condition of a constant 
density requires the mass scaling oc a^ . This corresponds 
to the initial mass ranging from w O.OOlAfg in the 3 AU 
case to K, 2OOM0 in the 200 AU case. With these values, 
all reference disks have about ten times the EKB density 
(10 EKBD). 

That all the belts share the same volume density of 
material is purely a matter of convention. Instead, we 
could choose them to have the same surface density or 
the same total mass. Given the scaling rules, as discussed 
in the text and Appendix [Aj none of these choices would 
have strong advantages or disadvantages. 

All five reference disks are listed in Table [TJ We evolved 
them with the collisional code, ACE, and stored all re- 
sults between the ages of 10 Myr and 10 Gyr at reason- 
able time steps. In what follows, we use self-explanatory 
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TABLE 1 
Description of reference disks 



Disk 
identifier 



Belt 
location [AU] 



Initial 
disk mass [Mg 



a range 

[AU] 



[AU] 



lOEKBD @ 3AU 


3 


0.001 


0.3 - 30 


0.5 - 20 


lOEKBD @ lOAU 


10 


0.03 


1- 100 


2-50 


lOEKBD @ 30AU 


30 


1 


3-300 


5-200 


lOEKBD @ lOOAU 


100 


30 


10 - 1000 


20 - 500 


lOEKBD @ 200AU 


200 


200 


20 - 2000 


30 - 1000 



identifiers like lOEKBD @ lOAU @ SOOMyr to refer to a 
particular disk of a particular age. 

Importantly, the same runs of the collisional code auto- 
matically provide the results for disks of any other initial 
density (o r mass). Thi s is po ssible due to the mass-time 
scaling of iLohne et al.l ()2008f ). which can be formulated 
as follows. Denote by M(Mq, t) the mass that a disk with 
initial mass Mq has at time t. Then, the mass of another 
disk with x times larger initial mass at time instant t/x 
is simply 

M{xMo, t/x) = xM{Mo, t), (4) 

For instance, the mass of the lEKBD @ lOAU @ lOGyr 
disk is one-tenth of the lOEKBD @ lOAU @ IGyr disk 
mass. Note that the same scaling applies to any other 
quantity directly proportional to the amount of disk ma- 
terial. In other words, M may equally stand for the mass 
of dust, its total cross section, thermal radiation flux, etc. 
See Appendix [Xl for additional explanations. 

4. RESULTS 
4.1. Size and spatial distributions of dust 

As noted above, the collisional code ACE uses masses 
and orbital elements of disk particles as phase space vari- 
ables. At any time instant, their phase space distribu- 
tion is transformed to usual mass/size and spatial dis- 
tributions. It is important to understand that mass/size 
distributions and spatial distributions cannot, generally, 
be decoupled from each other. Grains of different sizes 
have different radial distributions and conversely, the size 
distribution of material is different at different distances 
from the star. 

A typical size distribution of solids is shown in Fig. [3] 
for one of the disks, namely for lOEKBD @ 30AU @ 
lOOMyr. Different lines correspond to different distances 
from the primary. As expected, the size distribution is 
the broadest within the parent ring of planetesimals. Far- 
ther out, it only contains grains which are small enough 
to develop orbits with sufficiently large apocentric dis- 
tances due to radiation pressure. 

The spatial distribution of material in the same disk 
is shown in Fig. 3] Here, different lines refer to dif- 
ferent particle sizes. The ring of the biggest particles 
shown (100 /im), for which radiation pressure is negligi- 
ble, nearly coincides with the initial ring of planetesimals 
(semimajor axes: from 27 to 33 AU, eccentricities: from 
0.0 to 0.2, hence radial distances from 22 to 40 AU). The 
larger the particles, the more confined their rings. The 
rings are more extended outward with respect to the par- 
ent planetesimal ring than inward. 

Radial profiles of the normal geometrical optical depth 
for three reference disks (planetesimal rings at 10, 30, 
and 100 AU) are depicted in Fig. [5l Initially, the peak 
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Fig. 4. — Spatial distribution of grains with three character- 
istic radii for the lOEKBD @ 30AU @ lOOMyr disk. The ring of 
the biggest particles shown (100 /im, hatched) is the narrowest. Its 
radial extension is nearly the same as that of the initial planetes- 
imal ring; vertical "walls" are artifacts due to a discrete distance 
binning. 



optical depth of the disks is proportional to the distance 
of the parent ring, making the 100 AU disk ten times 
optically thicker than the 10 AU one. The subsequent 
collisional evolution of the disks depends on their initial 
mass and di stance from the star, as explained in detail in 
ILohne et al.l (|2008) and Appendix |^ Once a collisional 
steady state is reached (which is the case after 10 Myr for 
all three disks), the optical depth decays with time ap- 
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proximately as t^, where ^ « —0.3 ... — 0.4, i.e. roughly 
by one order of magnitude from 10 Myr to 10 Gyr. In 
a steady-state regime, the optical depth is proportional 
to r^+^-^i r^ r^-^ . This explains why, at any age between 



10 Myr and 10 Gyr, the 100 AU ring is 
cally thicker than the 10 AU one. 
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Fig. 5. — Radial profiles of the normal geometrical optical 
depth for three out of five basic runs (lOEKBD @ lOAU, solid 
lines; lOEKBD @ 30AU, dashed; lOEKBD @ lOOAU, dotted) at 
different ages. The thinner the line, the older the disk, as marked 
in the legend. The dashed-dotted lines are initial optical depths, 
artificially enhanced by a factor of ten for a better visibility. 



4.2. Dust temperatures 

Figure [H shows the dust temperatures as a function of 
two variables: grain distances from the star and their 
radii. In a parallel scale on the right, we show typical 
size distributions (cf. Fig. [3]). Similarly, under the tem- 
perature plot, typical radial profiles of the disk are drawn 
(cf. Fig. 11]). This enables a direct "read-out" of the typ- 
ical^ temperature in one or another disk. We find, for 
example, 130 K at 10 AU, 90 K at 30 AU, and 50 K at 
100 AU. 

These values are noticeably higher than the blackbody 
values of 88 K, 51 K, and 28 K, respectively. The reason 
for these big deviations and for the S-shaped isotherms in 
Fig. [6] is the astronomical sihcate's spectroscopic proper- 
ties with relatively high absorption at visible wavelengths 
and steeply decreasing absorption coefficient at longer 
wavelengths (see Fig. [2]). The cross-section dominat- 
ing astrosil grains are in a size range where the absorp- 
tion efficiency for visible and near-infrared wavelengths 
(around 1 fim) has already reached the blackbody value 
while emission is still rather inefficient. With the en- 
hancement of the emission efficiencies relative to the "sat- 
urated" absorption, temperatures drop drastically for 
somewhat larger grains. The larger the distance from the 
star (yielding lower average temperature and lower emis- 
sion efficiency), the wider the size range over which the 
temperature decreases, and the stronger the temperature 
difference between small and large grains. This explains 

^ "Typical" in the sense that it is the temperature of cross- 
section dominating grains in the densest part of the disk. 



why the S-shape of the isotherms gets more pronounced 
from the left to the right in Fig. [SI 

Further, we note that Mie resonances can increase the 
absorption/emission efficiencies even beyond unity for 
wavelengths somewhat longer than the grain size (see 
1 /im, 10 /im, and 100 /im curves in Fig.U). This explains 
the temperature maximum for grains of about 0.3 fim ra- 
dius ("resonance" with the stellar radiation maximum) 
and the minimum with temperatures even below the 
blackbody values for 10 to 50 /im grain radius ("reso- 
nance" with the blackbody emission peak). 

4.3. Spectral energy distributions 

We start with a single, "typical" SED for one of the 
disks. Such an SED for the lEKBD @ 30AU @ lOOMyr 
disk is shown in Fig. [7] with a thick solid line. It peaks at 
about 50 /im, which is consistent with the dust tempera- 
tures (Fig. [5]). The hump at « 10 /im is due to a classical 
silicate feature, as discussed below. 

For comparison, we have overplotted the SEDs calcu- 
lated for the same disk, but under different assumptions 
about the absorbing and emitting properties of grains: 
in a black-body approximation (grey line) and for amor- 
phous carbon (dashed line). Note that the difference 
applies only to the calculation of thermal emission. In 
other words, the dynamical modeling was still done by 
assuming the radiation pressure of astrosil and not of 
perfectly absorbing or carbon particles, but we assumed 
the grains to absorb and emit like a blackbody or carbon 
when calculating the thermal emission. There is a strik- 
ing difference between the curves, especially the black- 
body SED deviates from the others dramatically. The 
blackbody assumption leads to a strong increase of the 
total flux as well as to a shift of the maximum in the SED 
from 50 to 130 /im! In addition the excess drops towards 
longer wavelengths much slower than in the case of the 
astronomical silicate. In fact, it will never intersect the 
stellar photospheric flux. 

We now proceed with a set of SEDs for our grid of refer- 
ence disks. Some of them are shown in Fig. [8l The main 
features of these plots reveal no surprises. The absolute 
level of excess emission is higher for more massive disks, 
as well as for distant ones (which is just the consequence 
of the assumed "same-density" scaling, as described in 
Sect. 3.5, see also Fig. [5|). The amount of dust emission 
is roughly comparable with the photospheric emission for 
the mid-aged lEKBD @ 30AU disk. This is consistent 
with the known fact that a several Gyr-old EKB counter- 
part would only slightly enhance the photospheric emis- 
sion even at the "best" wavelengths. The position of the 
maximum emission ranges from « 30 /im for the lOAU 
disk to w 70 /im for the 100 AU disk. Note that black- 
body calculation would predict the emission to peak at 
longer wavelengths; beyond 100 /im for a 100 AU disk. 

Again, the hump seen in all SEDs slightly below 10 /im 
is due to a silicate feature in Qats] furthermore, some 
traces of the second feature at 20 /im are barely visible. 
This explanation is supported by Fig. [2] that shows the 
absorption efficiency feature in this spectral range for 
small particles. This becomes even more obvious by com- 
paring the contribution of the different grain size decades. 
For 0.1 to 1 /im particles the hump is more pronounced 
than for larger ones (see left panels in Fig. [5] below), as 
it is the case for the absorption efficiency. Further on. 
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Fig. 6. — The left upper plot shows the equilibrium temperature of dust particles as a function of their distance from the star (horizontal 
axis) and size (vertical axis). Contours are isotherms. The blackbody dust temperatures are given along the upper edge of the plot for 
comparison. In the right-hand plot the size distribution at the "central" distance of the systems (10 AU, solid; 30 AU, dashed; and 100 AU, 
dotted) at 100 Myr is given. The lovifcst left plot gives the normal optical depth for the same three disks as a function of distance to the 
star. An intersection of a horizontal straight line going through the maximum of the size distribution in a disk (right) with a vertical line 
through the peak of its radial profile (bottom) provides the typical dust temperature in that disk. 
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Fig. 7. — The emission from one and the same, lEKBD @ 30AU 
@ lOOMyr, disk, calculated under different assumptions about ab- 
sorbing and emitting properties of dust grains: blackbody, astrosil 
(our nominal case), and amorphous carbon particles. Thin solid 
line: photosphere of a G2V star. 



the 10 /im "excess" becomes less visible for most distant 
disks (from top to bottom panels in Fig. [5]), where the 
average temperatures arc lower, the maximum emission 
shifts to longer wavelengths, and therefore the Planck 
curve at A '^ 10 /um-20 ^m is steeper. 

The left panels in Fig. [9] illustrate relative contribu- 
tions of different-sized particles to the full SEDs. This 
is useful to get an idea which instrument is sensitive to 
which grain sizes. The blowout grains with radii less than 
0.4 /im make only modest contribution to the flux even 
at 10 /im. The mid-IR fluxes are always dominated by 
bound grains with 0.4 /im to 1 /im radii (for the 10 AU 
and 30 AU rings) or those with 0.4 /im to 10 /im (for the 
100 AU ring). In the far-IR, particles up to 100 /im in size 
play a role. The greatest effect on the sub-mm fluxes is 
that of 100 nm to 1 mm grains. 

The position of the different maxima in Fig. [9] can be 
understood by comparing the size decades to the dust 
temperature plot, Fig. [6] Particles of 0.1 /im to 1 /im 
are on the average a bit warmer than particles of 1 to 
10 /im. However, the size distribution shows that the sec- 
ond decade is dominated by particles only slightly larger 
than 1 /im, which are still nearly as warm as the parti- 
cles in the decade below. Thus, the maxima of the cor- 
responding SED contributions are shifted only slightly. 
It is the step to the next decade where the decrease of 
temperature becomes very obvious by a large shift of the 
maximum. From that size on, the maxima stay nearly at 
the same position (in fact the maxima are shifted again 
to smaller wavelengths) as the temperature changes only 
marginally. 

Similar to the contribution of the different size decades 
in the left panel, the right panels in Fig. [9] demonstrate 
the contribution of the different radial parts of the disk 
to the total SED. As expected, most of the flux comes 
from the medium distances as this is the location of the 
birth ring. The second largest contribution is made by 
the outer part of the ring. 

5. APPLICATION TO SELECTED DEBRIS DISKS 

5.1. Measured fluxes 

To test the plausibility of our models, we have selected 
several nearby sun-like stars known to possess debris 
dust. We used published datasets to search for stars 



Star Teff [K] log L./Lq D [pc] age [Myr] 



HD 377 
HD 70573 
HD 72905^ 
HD 107146 
HD 141943 



5852 " 
5841 " 
5831 " 
5859 " 
5805 " 



0.09 ") 
-0.23 "") 
-0.04 "' 
0.04 ") 
0.43 ") 



40 ") 
46 <") 
13.85 <*) 
29 °) 
67 ") 



32 "' 
100 ") 
420 ''^ 



100 



20 

32 ") 



References. 
IRhee at al.l iWOlh 
l|2008l ) 

Note. — ^A G1.5 star 



") [Hillenbrand et al." ("2001); *>' 
) IMoor et al., (,2006, 1. ,Trilling et al.l 



with (i) spectral classes most likely G2V (or very close), 
and (ii) unambiguous excesses probed in a wide range 
of wavelengths from near-IR to far-IR or sub-mm. The 
resulting five stars and their properties are listed in Ta- 
ble [21 a summary of observational data on them is given 
in Table [31 and the disk properties as derived in original 
papers are collected in Table [H The data include those 
from various surveys with IRAS, ISO, Spitzer, Keck II, 
and JCMT (Table [3]). The estimated ages of the systems 
range from 30 to 400 Myr (Table ^ and the fractional 
luminosities from ~10~^to--10~^ (Table [4]). The col- 
lected data points for our sample stars (photosphere -|- 
dust) are plotted in Fig. [TOl 

5.2. Observed excesses 

Symbols in Fig. 1 1 1 1 represent the observed excess emis- 
sion for our sample stars. In the cases where the pho- 
tospheric subtraction was done in the source papers, we 
just used the published data points. In the cases where 
only the total measured flux (star -I- dust) was given, we 
proceeded as follows. Three IRAC points (3.6, 4.5, and 
8.0 /im) were fitted by an appropriate NextGen model 
(jHauschildt et al.lll999f ). and the resulting photosphcric 
spectrum was subtracted from the fluxes measured at 
longer wavelengths. As far as the data quality is con- 
cerned, the best case is clearly HD 107146, where the 
data points cover a broad range between 10 /im and 
1 mm. In other cases, the longest wavelengths probed 
lay at 70-160 /im. As a result, it is sometimes unclear 
where exactly the excess peaks. This is exemplified by 
HD 70573 where the 160 /im point has a huge error bar. 

Yet before any comparison with the modeled SEDs, 
the resulting points in Fig. llll allow several quick conclu- 
sions. Notwithstanding the paucity of long-wavelength 
data just discussed, in all five systems the excess seems 
to peak at or slightly beyond 100 /im, suggesting a "cold 
EKE" as a source of dust. Additionally, in all systems 
except for HD 377, a warm emission at A < 20 /im seems 
to be present, implying a closer- in "asteroid belt". 

5.3. Comparison of measured and modeled SEDs 

We now proceed with a comparison between the ob- 
served dust emission and the modeled emission. We 
stress that our goal here is not to provide the best fit to 
the observations possible with our approach, but rather 
to demonstrate that a set of reference disks modeled in 
the previous sections can be used to make rough prelim- 
inary conclusions about the planetesimal families. 
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TABLE 3 

Observational data for the five G2 stars and their disks 



Star 



Instrument, A (^in) 



Reference 



Notes 



HD 377 



HD 70573 



ffiAC 3.6/4.5/8.0 
IRAS 13/33 
IRAS 60 
MIPS 24/70/160 



J«AC' 3.6/4.5/8.0 
IRS 13/33 
MIPS 24/70/160 



Hillenbrand et al. (2008) 
Hillenbrand et al. (2008) 
Moor et al. (2006) 
Hillenbrand et al. (2008) 



Hillenbrand et al. (2008) 
Hillenbrand et al. (2008) 
Hillenbrand et al. (2008) 



A planet host star 
(Setiawan et al. 2007) 



ffiAC 3.6/4.5/8.0 
IRS 13/33 
IRAS 12/25 
JSOPHOT 60/90 
MIPS 24 
MIPS 70 



Hillenbrand et al. ( 201)87 
Beichman et al. (2 006) 
Spangler et al. (20^7 
Spangler et al. (20W) 
Brvden et al. (2006) 
Hillenbrand et al. (^008) 



HD 72905 



HD 107146 /iJAC 3.6/4.5/8.0 Hillenbrand et al. (2008) 

LWS 11.7/17.8 Williams et al. (2004) 

IRS 13/33 Hillenbrand et al. (2008) 

IRAS 60/100 Moor et al. (2006) 

MIPS 24/70 Hillenbrand et al. ( 2008) 

SCUBA 450/850 Williams et al. (2004f~ 



HD 141943 ffiylC 3.6/4.5/8.0 Hillenbrand et al. (2008) 

IRS 13/33 Hillenbrand et al. (2008) 

MIPS 24/70 Hillenbrand et al. (2008) 



Resolved in V and I 
bands (Ardila et al. 2004) . 
at 350 and 450 /xm 
(Williams et al. 2004') . 
and at 3 mm 
(Carpenter et al. 2005) 



TABLE 4 
Previously derived disk properties 



Star 



[K] 



Rd^st [AU] 



M. 



[M(&] 



-^dust/-^* 



HD 377 



58 "'•! 



23 "'-^ 



3.98 X 10-* ")• 



3.98 X 10 



-4 a), 11 



35 ""''S 



(4.0 ±0.3) X 10 



-4 /),12 



4 a), 11 



HD 70573 



41 °-> 



2.0 X 10" 



1.0 X 10^ 



7 a), 5 

12.2- 15.9 ''''3 
6.2 9)'5 



1.58 X lO"'^ "")•* 
3.3 X lO-*^ ')'3 



-5 a), 11 
-5 i)),13 
lQ-5 9), 14 
-5 e),15 
-4 9), 16 



HD 72905 103 "^'^ 

63 - 67 *''^ 
123 9)'2 



2.0 X 10" 
2.9 X 10" 
(0.6- 1.5) X 
1.6 X 10" 
2.8 X 10 



JTX 



30 °-'>'^ 
13.6- > 200 ")■'' 

29 d),5 
> 31 - 150 '•'■^ 



1.26 X 10-3 ")■* 4.94 X 10"* '"' 



11 



HD 107146 



52 



55 '')'2 
51 '').* 



3.2 X lO-'^ "'''^ 

8.99 X 10-2 '*)'io 

0.1 '')'4 



J7X 



18 "J'S 
8.6-40 ")■'' 



7.94 X 10-^ "'■** 



(9.2 ±0.9) X 
9.5 X 10- 
1.2 X 10 



10-4 /),12 
4 d),12 
-3 h) 



1.58 X 10-* "^'^^ 



HD 141943 



85 



References. — - a) [H illenbrand et al. (2008), b) Beiehman et al. (2006), c) 
[Carpenter et al.l 1120051). d) [Rhee et al., (i2007i), e) iBrvden et al.i (i2006i ). f ) iMoor et al.i (i2006i ). 
g) ISpangler et al.l II2001D . h) IWilliams et al.l EOM) 



Note. — ^ Color temperature (33 - 70/im) from blackbody SED fitting. ^ From SED fitting 
using a single temperature blackbody. From SED fitt ing using 10/im silicat e grains with a 

temperature profile following a power law (favored model in lBeichman et al.l < I2006I) ). From single 
temperature SED fitting using a modified blackbody and a mass absorption coefficient K85o — 
1.7 cm /g. Derived from T^ust assuming blackbody (lower limit). Extended ring derived 

from blackbody SED fitting assuming a constant surface density. Inner border derived from SED 
fitting, outer border taken from resolved image. Derived from fractional luminosity for an average 
grain size of < a >— 10/^m and a density of p — 2.5g/cni . Derived for T^ugt — 40 K using 
a frequency dependent mass absorption coefficient. Derived from submillimetcr observations 

using a dust opacity of 1.7 cm /g at 850 /im. Derived from T^ust arid i?.duat using Stefan- 

Boltzmann relation. Ldust/L* — Lin/ Lhol- .^dust obtained by integrating IRS spectrum 

(10 — 34 fim) after extrapolation to 70 fim. .^dust is derived from the SED fitting and L^ 

is obtained by integrating the corresponding Kurucz model. Minimum value, derived from 

the 70 fim measurement. L« is the stellar bolometric luminosity and Ldust is the sum of the 

luminosities in each {IRAS) wavelength band with a correction (for longer wavelengths). 
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Fig. 8. — Spectral energy distributions of disks stemming from planetesimal rings with different masses at different locations and at 
different time steps. To obtain the absolute values of fluxes, a distance of fO pc was assumed. Left: reference disks (lOEKBD), right: 
less massive disks (lEKBD). The results for the latter have been obtained with the aid of Eq. l|4]l. From top to bottom: the SEDs of the 
simulated planetesimal rings at 10, 30 and 100 AU. In each panel, lines of decreasing thickness correspond to the ages of 10 Myr, 100 Myr, 
1 Gyr, and 10 Gyr. Note that the evolution of the lEKBD @ lOOAU disk at the beginning is very slow, so that the SEDs at 10 and 100 Myr 
are indistinguishable. Vertical lines indicate centers of observational bands of several instruments (in fira): Spitzer MIPS (24, 70, 160), 
Herschel PACS (100, 160), Sofia HAWC (200), CSO Share (350), JCMT SCUBA/SCUBA 2 (450, 850), MPIfR, IRAM (1300). A thin line 
from top left to bottom right is the stellar photosphere. 
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Fig. 9. — Contribution of individual grain size decades (shown with different hnestyles in the left panel) and individual radial annuli of 
the disks (different Hnestyles, right) to the SED. As the grain blowout radius is ss 0.4 /xm, see Fig. [T] in the left panels we split the lowest 
size decade into blowout grains with s € [0.1 /.im, 0.4^m] and bound ones with s G [0.4/im, 1.0/.tm]. Panels from top to bottom correspond 
to planetesimal rings at 10, 30 and 100 AU. The initial density of all disks is lOEKBD and their age is 1 Gyr. 



To make such a comparison, we employ the following 
procedure: 

1. For each star, we first look whether only cold or 
cold + warm excess emission is present. In the former 
case (HD 377), we fit the data points with a single "cold" 
reference disk. In the latter case (all other systems), we 
invoke a two-component model: a close-in 3 AU disk and 
an appropriate "cold" disk. 

2. The location of the "cold" planetesimal belt is cho- 
sen according to the peak wavelength of the measured 
excess: 100 AU (HD 72905 and HD 141943) or 200 AU 
(HD 377, HD 70573, and HD 107146). 

3. We then scale each of the two reference SEDs, 
"warm" and "cold" (or only one for HD 377) vertically to 
come to the observed absolute flux. Physically, it neces- 
sitates a change in the initial disk mass. However, it is 



not sufficient to change the initial disk mass by the ratio 
of the observed flux and the flux from a reference disk. 
The reason is that a change in the initial mass also alters 
the rate of the coUisional evolution, whereas we need the 
"right" flux at a fixed time instant, namely the actual 
age of the system (Tab. [2|). Therefore, to find the mass 
modification factor we apply scaling rules, as explained 
in Appendix [Al Specifically, we solve Eq. (jA8|l . In the 
systems that reveal both warm and cold emission, this is 
done separately for the inner and outer disk. 

The results presented in Fig. [Tl] with lines show that 
the modeled SEDs can, generally, reproduce the data 
points within their error bars. Again, the judgment 
should take into account the fact that we are just us- 
ing one or two pre-generated SEDs for a rather coarse 
grid of reference disks. Much better fits would certainly 
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X[|j,m] 

Fig. 10. — Observational data for five selected G2V stars. Note that all fluxes have been scaled to the same standard distance of 
lOpc. Symbols in the left-hand, grey-shaded part of each panel (A < 10 fim) are IRAC observations. They are used to find an appropriate 
Hauschildt model to the photosphere (thin solid line), assuming that no excess is already present in the near infrared. Vertical error bars 
are la observational uncertainties, taken from the source papers. Horizontal bars indicate the band width of the respective detector. 
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be possible if we allowed a more exact positioning of par- 
ent belts and let additional model parameters vary. Dust 
opacities, initial distributions of planetesimals' sizes and 
orbital elements, as well as their mechanical properties 
that were fixed in modeling of the coUisional outcomes 
would all be at our disposal for this purpose. Further, 
more than two-component planetesimal belts could be as- 
trophysically relevant as well, as is the case in our solar 
system (asteroid belt, different cometary families, vari- 
ous populations in the EKB). 

We now come to the interpretation of the fitting re- 
sults, trying to recover the properties of dust-producing 
planetesimal belts. Table [5] lists them for all systems. 
The most important information is the deduced mass 
and location of the belts. 

5.4. Results for hot dust 

As far as the hot dust components in four out of five 
systems are concerned, our results show that these can 
be explained by "massive asteroid belts" with roughly 
the lunar mass in bodies up to ~ 100 km in size, lo- 
cated at 3 AU, with a width of ~ 1 AU. However, the 
quoted distance of inner components — 3 AU — is only 
due to the fact that this is the smallest disk in our grid. 
This distance can only be considered as an upper limit: 
the SEDs seem perfectly compatible with disks as far 
in as 0.3 AU, as suggested for the case of HD 72905 
(|Wvatt et al.l[200l . 

What is more, even the very fact that hot excess is real 
can sometimes be questioned, since it can be mimicked by 
photospheric emission slightly larger than the assumed 
values. Indeed, the excess for HD 70573 and HD 72905 
at wavelengths around and below 25 /im does not ex- 
ceed 10%, which is comparable with the average cali- 
bration uncertainty and therefore h as to be considere d 
marginal (Bryden et al. 2006; Hille nbrand et all l2008l ). 
Only in the case of HD 72905, the Spitzer/IRS detec- 
tion of the 10 /Ltm emission from hot silicates provides 
an independent c onfirmation that the hot excess is real 
(jBeichma n et all [2006). However, the HD 72905 plot in 
Fig. [Tl] makes it obvious that some problems occurred 
in terms of the photosphere fitting. All data points that 
we obtained by subtracting the IRAC photospheric fluxes 
(squares) systematically lie above the data points where a 
photosphere from the literature was subtracted (circles). 
The origin of the difference is unclear; on any account, 
the problem cannot be mitigated by the assumption that 
an excess is already present at IRAC wavelengths, since 
this would shift the squares further upwards. Consider- 
ing the circles to be more trustworthy, the shape of the 
SED to fit changes. Then a closer-in disk at ^ 0.3 AU 
could better reproduce the fluxes in the near and mid 
infrared, while the outer ring would have to be shifted to 
a distance somewhat larger than 100 AU in order not to 
surpass the measured flux at 33 /um. A problem would 
arise with the inner disk: at ~ 0.3 AU, the coUisional evo- 
lution is so rapid that an unrealistically large initial belt 
mass would be ne cessary. Similar arguments have led 
IWvatt et all ()2007f ) to a conclusion that HD 72905 must 
be a system at a transient phase rather than a system 
collisionally evolving in a steady state. 

Still, treating the derived sizes and masses of the inner 
disks as upper limits yields physical implications. Be- 
cause the coUisional evolution close to the star is rapid. 



such belts must have lost up to two-thirds of their ini- 
tial mass before they have reached their present age (cf. 
initial and current mass in Table [S]). In the case of 
HD 70573, t he known giant plan et with a = 1.76 AU 
and e = 0.4 (jSetiawan et al.l '2007l does not seem to ex- 
clude the existence of a dynamically stable planetesimal 
belt either inside ~ 1 AU or outside ^ 3 AU. 

5.5. Results for cold dust 

The estimated parameters of the outer components of 
the disks suggest "massive and large Kuiper belts" . The 
radii of the outer rings are larger than the radii derived 
in previous studies (cf. Table [4] and Table [5]). This traces 
back to our using astrosilicate instead of blackbody when 
calculating the dust emission, so that the same dust tem- 
peratures are attained at larger distances (see Fig. [7]) . 

Since one disk in our sample, that of HD 107146, 
has been resolved, it is natural to compare our de- 
rived disk radius with the one obtained from the images. 
' Williams et al.l (2004) report an outer border of the sys- 
tem o f 150 AU bas e d on submillimeter images. In con- 
trast, lArdila et all (|2004D detected an 85 AU- wide ring 
peaking in density at about 130 AU. This is compa- 
rable to, although somewhat smaller than, our 200 AU 
radius. However, moving the outer ring to smaller dis- 
tances would increase the fluxes in the mid infrared 
where the SED already surpasses the observations and 
the other way round in the sub-mm region. The re- 
sulting deficiency of sub-mm fluxes, though, could be 
du e to roughness of Mie c alculations. As pointed out 
bv IStognienko et all ()1995l ). an assumption of homoge- 
neous particles typically leads to underestimation of the 
amount of thermal radiation in the sub-mm region. 

Large belt radii imply large masses. Dust masses de- 
rived here are by two orders of magnitude larger than 
previous estimates (cf. Table [5] and Table |4|). The total 
masses of the belts we derive range from several to sev- 
eral tens earth masses, to be compared with ~ O.IM^ 
in the present-day EKB (although there is no u nanim- 
ity on that point — cf. IStern fc Colwel]||1997t ). Note 
that, as the coUisional evolution at 100-200 AU is quite 
slow, whereas the oldest system in our sample is only 
420 Myr old, the difference between the initial disk mass 
and the current disk mass is negligible. Assuming several 
times the minimum mass solar nebula with a standard 
surface density of s olids E ^ 50gcm^(r/l AU)^'^/^ (e.g. 
iHavashi et al.l [19851 ). the mass of solids in the EKB re- 
gion would be a few t ens M^; and current models (e.g. 
iKenvon &: Luul [r999bD successfully accumulate 100 km- 
sized EKB objects in tens of Myr. However, it is ques- 
tionable whether the assumed radial surface density pro- 
file could extend much farther out from the star. As a 
result, it is difficult to say, whether a progenitor disk 
could contain enough solids as far as at 200 AU from the 
star to form a belt of 30 to SOMq. 

However, such questions may be somewhat premature. 
On the observational side, more data are needed, espe- 
cially at longer wavelengths; for instance, the anticipated 
Herschel data (PACS at 100/160 /xm and SPIRE at 250 
to 500 /im) would help a lot. On the modeling side, a 
more systematic study is needed to clarify, how strongly 
various assumptions of the current model (especially the 
coUisional outcome prescription and the material choices) 
may affect the calculated size distributions of dust, the 
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TABLE 5 
Disk properties derived in this study 



Star 


Component 


Mdisk [A%] 1) 


Rbcit [AU] 2) 


A/du.t [Me] 3) 


Tdu.t [K] 4) 


HD 377 


Outer 


(32) 


32 


200 


3.1 X 10-2 


40 


HD 70573 


Inner 
Outer 


(0.0063) 
(2.6) 


0.0046 
2.5 


3 

200 


1.4 X 10-' 
2.0 X 10-3 


200 
40 


HD 72905 


Inner 
Outer 


(0.054) 
(0.23) 


0.019 
0.23 


3 
100 


3.4 X lO-'' 
2.1 X 10-* 


200 
50 


HD 107146 


Inner 
Outer 


(0.039) 

(47) 


0.023 

47 


3 

200 


4.9 X 10- "■ 
4.8 X 10-2 


200 

40 


HD 141943 


Inner 
Outer 


(0.039) 
(6.1) 


0.027 
6.1 


3 
100 


8.0 X 10-' 
5.5 X 10-3 


200 
50 



Note. — ^' Initial mass (in parentheses) and the current mass of the whole planetesimal 
disk (bodies up to 100 km in radius). 

2) Location of the parent planetesimal belt. 

3) Current mass of "visible" dust (grains up to 1mm in radius). 

*) Temperature of cross-section dominating astrosil grains at the location of the parent plan- 
etesimal belt, see explanation at Fig. \E\ 



dust grain temperatures, and the amount of their ther- 
mal emission. 

At this point, we can only state that in the five sys- 
tems analyzed (with a possible exception of HD 72905) 
and with the caveat that available data are quite scarce, 
the observations are not incompatible with a stan- 
dard steady-state scenario of collisional evolution and 
dust production. Of course, other possibilities, such 
as major collisional br eakups (jKenvon &: BromlevI 120051 : 
iGrigorieva et al.1 I2007D or events similar to the Late 
Heavy Bombardment (as suggested, for instance for HD 
72905, IWvatt et all l2007t l cannot be ruled out for the 
inner disks. 

6. SUMMARY 

Debris disks around main-sequence stars may serve as 
tracers of planetesimal populations that have accumu- 
lated at earlier, protoplanetary and transitional, phases 
of systems' evolution, and have not been used up to form 
planets. However, observations of debris disks are only 
sensitive to the lowest end of the size distribution. Using 
dynamical and collisional models of debris disks is the 
only way to "climb up" the ladder of the collisional cas- 
cade, past the ubiquitous /im-sized grains towards parent 
bodies and towards the main mass reservoir of the disks. 

The main idea of this paper has been to take a grid of 
planetesimal families (with different initial masses, dis- 
tances from a central star etc.), to collisionally "gener- 
ate" debris disks from these families and evolve them 
with the aid of an elaborated collisional code, and finally, 
to calculate SEDs for these disks. A comparison/fit of the 
observed SEDs with the pre-generated SEDs is meant to 
allow quick conclusions about the properties of the plan- 
etesimal belt(s) that maintain one or another observed 
disk. 

Our specific results are as follows: 

1. We have produced five reference disks around a G2V 
star from planetesimal belts at 3, 10, 30, 100, and 200 AU 
with 10 times the EKB mass density and evolved them 
for 10 Gyr. With an appropriate scahng rule (Eg. lAip . 
we can translate these results to an arbitrary initial disk 
mass and any age between 10 Myr and 10 Gyr. Thus, ef- 
fectively we have a three-parametric set of reference disks 
(initial mass, location of planetesimal belt, age). For all 



the disks, we have generated SEDs, assuming astrosili- 
cate (with tests made also for blackbody and amorphous 
carbon). 

2. We have selected five G2V stars with good data 
(IRAS, ISO/ISOPHOT, Spitzer/IRAC, /IRS, /MIPS, 
Keck II/LWS, and JCMT/SCUBA) and tested our grid 
against these data. For all five systems, we have repro- 
duced the data points within the error bars with a lin- 
ear combination of two disks from the grid (an "asteroid 
belt" at 3AU and an outer "Kuiper belt"). This auto- 
matically gives us the desired estimates of planetesimals 
(location, total mass etc.). 

3. A comparison of the observational data on the five 
stars with the grid of models leads us to a conclusion 
that the cold emission (with a maximum at the far-IR) 
is compatible with "large Kuiper belts" , with masses in 
the range 3-50 earth masses and radii of 100-200 AU. 
These large sizes trace back to the facts that the colli- 
sional model predicts the observed emission to stem from 
micron-sized dust grains, whose temperatures are well in 
excess of a blackbody temperature a t a given distance 
from the star (as discussed, e.g., in Hillenb rand et aLl 
l2008f l. This conclusion is rather robust against variation 
in parameters of collisional and thermal emission models, 
and is roughly consistent with disk radii revealed in scat- 
tered light images (e.g. HD 107146). Still, quantitative 
conclusions about the mass and location of the planetes- 
imal belts would significantly depend on (i) the adopted 
model of collision outcomes (which, in turn, depend on 
the dynamical excitation of the belts, i.e. on orbital ec- 
centricities and inclinations of planetesimals) and (ii) the 
assumed grains' absorption and emission efficiencies. For 
example, a less efficient cratering (retaining more grains 
with radii ~ 10 /zm in the disk) and/or more "transpar- 
ent" materials (making dust grains of the same sizes at 
the same locations colder) would result in "shifting" the 
parent belts closer to the star. 

In future, we plan to extend this study in two direc- 
tions. First, we will investigate more systematically the 
influence of the dust composition by trying relevant ma- 
terials with available optical data rather than astrosil- 
icate; this should be done consistently in the dynami- 
cal/collisional and thermal emission models. Second, it 
is planned to extend this study to stars with a range 
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Fig. 11. — Observed (s ym bols) and modeled (lines) exce ss emission, scaled to the distance of lOpc. The wavelength range matches 
the unshaded part of Fig. 1101 Here, in contrast to Fig. 1101 symbols represent the excess emission. Squares mark the cases where the 
scaled NextGen model shown on that figure was used to subtract the photosphere. Circles indicate that for these observations the stellar 
emission was subtracted using photospheric fluxes as given in the respective papers. Dashed lines: two "underlying" SEDs of reference 
disks (unsealed, i.e. with lOEKBD), one for "cold" excess and one for "warm" excess (except for HD 377 where only cold component is 
observed). SoUd fine: a Hnear combination of two scaled reference SEDs that provides a reasonable fit to the data points (except for HD 377 
where a single scaled reference SED is sufflcient). 
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of spectral classes. This will result in a catalog of disk 
colors that should be helpful for interpretation of data 
expected to come, most notably from the Herschel Space 
Observatory. 

This work has been particularly motivated by the Her- 
schel Open Time Key Program "DUNES" (DUst disks 
around NEarby Stars, PI: C.Eiroa) and we wish to 
thank many colleagues involved in DUNES (in particular, 



Jean-Charles Augereau, Jens Rodmann, and Philippe 
Thebault) for encouragement and numerous discussions. 
A speedy and constructive review of an anonymous re- 
viewer helped to improve the paper. This research 
has been funded by the Deutsche Forschungsgemein- 
schaft (DFG), projects Kr 2164/5-1 and Mu 1164/6-1, by 
the Deutscher Akademischer Austauschdienst (DAAD), 
project D/0707543, and by the International Space Sci- 
ence Institute (Bern). 



APPENDIX 
SCALING RULES 



1. Dependence of evolution on initial disk mass. Consider a disk with initial mass M{t = 0) = Mq at a distance r 
from the star with age t. Denote by F(Mq, r, t) any quantity directly proportional to the amount of disk material in 
any size regime, from dust grains to planetesimals. I n other words, F m ay equally stand for the total disk mass, the 
mass of dust, its total cross section, etc. As found by lLohne et al.l ()2008[ ). there is a scaling rule: 



F{xMo,r, t) = xF{MQ,r, xt), 



(Al) 



valid for any factor x > 0. This scaling is an exact property of every disk of particles, provided these are produced, 
modified and lost in binary collisions and not in any other physical processes. 

2. Dependence of evoluti on on distance . Another sc aling rule is the dependence of the evolution timescale on the 
distance from the star (Wvatt et aT]|2007l: lL5hne et al.i 20081 . Then 



F{Mo,xr,t)^F{Mo,r,t-^-n 



(A2) 



Unlike Eq. (JA1|1 . this scaling is approximate. 



3. Dust mass as a function of time. Finally, the third scaling rule found in iLohne et al.l (|2008[ ) is the power-law 
decay of the dust mass 

F(Mo, r, xt) » x-^F{Mo, r, i), (A3) 



where ^ w 0.3 . . . 0.4 (Fig. [T^ . This scaling is also approximate and, unlike Eq. (lAip and Eq. (|A2p . only applies to 
every quantity directly proportional to the amount of dust. In this context, "dust" refers to all objects in the strength 
rather than gravity regime, implying radii less than about 100 meters. The scaling is sufficiently accurate for disks 
that are much older than the coUisional lifetime of these lOOm-sized bodies. This is also seen in Fig. [121 while for the 
3 AU disk the power law (|A3p sets in after <C 1 Myr, the 200 AU disk needs ~ 100 Myr to reach this regime. 
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Fig. 12. — The time evolution of dust mass {s < 1mm) for our five reference disks (thick fines). For comparison, power laws t ^ with 
^ = 0.3 and ^ = 0.4 are shown with thin dashed lines. 



Note that the "pre-steady-state" phase of collisional evolution may actually require a more sophisticated treatment. 
Our runs assume initially a power-law size distribution of planetesimals, and an instantaneous start of the collisional 
cascade at i = 0. In reality, an initial size distribution is set up by the accretion history of planetesimals and will surely 
deviate from a single power law. Moreover, at a certain phase cratering and destruction of objects may increasingly 
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come into play simultaneously with ceasing, yet ongoing accretion; the efficienc ies and timescales of these processes will 
be di f ferent for different size ranges and different spatial locales in the disk (e.g. lDavis fc Farinellal[l997t iKenvon fc Luul 
fl998lfl999a llbll. 

The usefulness of these scaling rules can be illustrated with the following examples. 



Example 1. Assuming now F to be the total amount of dust, from Eqs. (|A1|1 - (|A3P one finds 

F{xM^, yr, t) « x^-^y^-^^^FiM^, r, t). (A4) 

Our choice of reference disks (different distances, but the same volume density) implies x = y^. The normal optical 
optical depth r ex F/r'^ scales as 

T{yH'Io,yr,t) « 2/i+i-3«T(Mo,r,t). (A5) 

Therefore, once a steady-state is reached (^ « 0.3. ..0.4), a y times more distant planetesimal belt gives rise to a y^+^-^i 
times optically thicker disk. This explains, in particular, why in Fig. [4] any 100 AU ring is « 30 times optically thicker 
than the co-eval 10 AU one. 



Example 2. Since the distance r in Eqs. (JA1[) and (jA3[) is kept fixed, F in these equations can also denote the 
radiation flux, emitted by a disk at a certain wavelength. Let Fo{t) be the observed fiux from a disk of age t. Imagine 
a model of a disk of the same age with an initial mass Mq predicts a flux F{Mq, t) which is by a factor A lower than 
the observed one: 

Fo{t)^AF(Mo,r,t). (A6) 

Our goal is to find the "right" initial mass, i.e. a factor B such that 



^o(i)=F(BAfo,r,t). 
With the aid of Eq. (jAip . this can be rewritten as 

Fo{t) = BF{Mo,r,Bt). 



Fo{t) « BF{Mo,r,t)B-^ = B^-^F(M^,r,t), 



(A7) 

(A8) 
(A9) 

(AID) 

For instance, a 10 times higher flux at a certain age requires a 27-46 times larger initial disk mass if ^ = 0.3 . . . 0.4. 

Although this rule is convenient for quick estimates, it should be used with caution. As described above, the value 
of ^ at the beginning of coUisional evolution (which lasts up to 100 Myr for the 200 AU belt) can be much smaller — 
close to zero or even negative — than the "normal" ^ — 0.3 . . . 0.4. For this reason, we prefer to use only the first 
scaling rule, Eq. (|A1[) . Therefore, instead of applying Eq. (|A10p . we find B by solving Eq. (jA8[) numerically with a 
simple iterative routine. It is this way Fig. 1111 was constructed. 



Eq. ()A3P gives now 



whence 



i?«Ai/(i-«). 
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